Tunneling of ultracold atoms in time-independent potentials 
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We present theoretical as well as experimental results on resonantly enhanced quantum tunneling 
of Bose-Einstein condensates in optical lattices both in the linear case of single particle dynamics and 
in the presence of atom-atom interactions. Our results demonstrate the usefulness of condensates 
in optical lattices for the dynamical control of tunneling and for simulating Hamiltonians originally 
used for describing solid state phenomena. 
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Tunneling as a quantum mechanical effect takes place in a classically forbidden region between two regions of clas- 
sically allowed motion. While the term "dynamical tunneling" typically refers to tunneling of quantum states across 
c/j'' dynamical barriers in classical phase space the original problem simply intended tunneling across a potential 
barrier. Both types of tunneling are addressed in this chapter, with major focus on situations in which external 
forces make the studied systems intrinsically time-dependent and allow for a dynamical control of tunneling through 
potential barriers or across band gaps which are dynamically explored by the system. 

A standard example of tunneling across static barriers is the motion in a double well potential. The two potential 
^2 i wells are separated by a potential barrier which is impenetrable for a low-energy classical particle. The quantum 
qh. mechanical solution shows that the wave packet initially localized in one of the wells performs oscillations between 
the two classically allowed region. Tunneling takes place between two levels nearly degenerate in energy, and in most 
Ci , cases the investigated tunneling takes place between the lowest energy states, for instance of a double well. However, 
^ in a potential configuration as the asymmetric double well shown in Fig. [2a), an energy matching between a ground 
, state on one side and an excited state on the other side leads to a tunneling between those states resonantly enhanced 
by the energy matching. In the resonantly enhanced tunneling (RET) the probability for the quantum tunneling of a 
particle between two potential wells is increased when the energies of the initial and final states of the process coincide. 
^ I In the one-dimensional double potential barrier of Fig. [IJb), the narrow central potential well has weakly-quantized 
—'I (or quasi-stationary) bound states, of which the energies are denoted by Ei and E2 in Fig. [TJ If the energy E of 
electrons incident on the barrier coincides with these energies, the electrons may tunnel through both barriers without 
any attenuation. The transmission coefficient reaches unity at the electron energy E — Ei 01 E — E2. It is interesting 
that while the transmission coefficient of a potential barrier is always lower than one, two barriers in a row can be 
' completely transparent for certain energies of the incident particle. 
[ In the early 1970s, Tsu, Esaki, and Chang computed the two terminal current-voltage characteristics of a fi- 

■ nite superlattice, and predicted that RET to be observed not only in the transmission coefficient but also in the 
current- voltage characteristic 0, Resonant tunneling also occurs in potential profiles with more than two barriers. 

, Technical advances led to the observation of negative differential conductance at terahertz frequencies and triggered a 
T-H ■ considerable research effort to study tunneling through multi-barrier structures. Owing to the fundamental nature of 
, this effect and the practical interest [1] , in the last few years much progress has been made in constructing solid state 
^ ■ systems such as superlattices [l-0|i quantum wells [1], and waveguide arrays which enable the controlled obser- 
vation and application of RET. The potential profiles required for resonant tunneling and realized in semiconductor 
system using heterojunctions allowed the manufacture of resonant tunneling diodes. These devices have important 

■ applications such as in high-frequency signal generation and multi- valued data storage, as reviewed in |l(^ . 

■ ■ In the last decade, the experimental techniques used in atom and quantum optics have made it possible to control 
the external and internal degrees of freedoms of ultracold atoms with a very high degree of precision. Thus, ultracold 
bosons or fermions loaded into the periodic optical potential created by interfering laser beams (double-well, lattices 
and superlattices) are optimal realizations of quantum mechanical processes and phenomena proposed and studied in 
other contexts of solid-state physics. Ultracold atoms and Bose-Einstein condensates, for instance, have been used 
to simulate phenomena such as Bloch oscillations in tilted periodic potentials [llHlsj and to study quantum phase 
transitions driven by atom-atom interactions [l6j . 

RET- like effects have been observed in a number of experiments to date. In ref. (iTj . resonant tunneling was 
observed for cold atoms trapped by an optical lattice when an applied magnetic field produced a Zeeman splitting 
of the energy levels. Resonant tunneling has been observed in a Mott insulator within an optical lattice, where a 
finite amount of energy given by the on-site interaction energy is required to create a particle-hole excitation [Tsj . 
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FIG. 1: (a) schematic representation of the energy levels within an asymmetric double well. The unperturbed energies within 
the left and right wells are indicated by the continuous lines. Because of the resonant tunneling between the ground state in 
the left well and the first excited one in the right well, the asymmetric and antisymmetric states have energy indicated by the 
dashed lines, (b) schematic band diagram of a resonant-tunneling diode structure under a voltage bias between the incoming 
(left) and outgoing (right) regions. 



Tunneling of the atoms is therefore suppressed. If the lattice potential is tilted by application of a potential gradient, 
RET is allowed whenever the energy difference between neighboring lattice sites due to the potential gradient matches 
the on-site interaction energy. This RET control in a Mott insulator allowed Foiling et al. to observe a second-order 
coherence, i.e. a two-atom RET (19j . 

Most of the quantum transport phenomena investigated with Bose-Einstein condensates within periodic optical lat- 
tices focused on the atomic motion in the ground state band of the periodic lattice. Only a few experiments examined 
the quantum transport associated with interband transitions "vertical" in the energy space. Interband transitions 
were induced by additional electromagnetic fields, as in the case of the spectroscopy of Wannier-Stark levels [20|, or by 
quantum tunneling between the bands. Tunneling between otherwise uncoupled energy bands occurs when the bands 
are coupled by an additional force, which can be a static Stark force (tilting the otherwise periodic lattice) or 
also by strong atom-atom interactions as observed for fermions in f2l'| and discussed for bosons in [25] . The quantum 
tunneling between the ground and the first excited band is particularly pronounced in the presence of degeneracies of 
the single- well energy levels within the optical lattice leading to RET. In [2^ [2^ such a type of RET was investigated 
for a Bose-Einstein condensate in a one-dimensional optical lattice, which allows for a high level of control on the 
potential depth and the lattice tilt. Those experimental investigations concentrated on the regime of parameters for 
which the tilting force - at RET conditions equal to the energy difference between neighboring wells - dominated the 
dynamics of the condensate. The RET tunneling of the ground band and the first two excited energy bands were 
measured in a wide range of experimental conditions. In addition the RET process is modified by the atom-atom 
interactions, bringing new physics to the quantum tunneling. 

This chapter is organized as follows. Section II sets the stage discussing optical lattices and giving the necessary 
background. While Section III reports on RET in closed two and three well systems. Section IV focuses on our main 
subject, the control of tunneling by RET in open quantum systems. This Section reports on our experimental data 
in the linear tunneling regime, i.e. in the absence of atom-atom interactions, as well as on interaction induced effects. 
In Section V a model for many-body tunneling is introduced, before we summarize the recent advances concerning 
RET in Section Ell 
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FIG. 2: (a) In an optical lattice without additional external forces, the ground-state levels are resonantly coupled, leading to a 
tunneling energy J. (b) When a linear potential is applied, e.g. by applying a force F, the levels are shifted out of resonance 
and tunneling is suppressed (Wannier-Stark localization), (c) If an additional potential energy oscillating at an appropriate 
frequency u) is applied, the levels can again be coupled through photons of energy hw and tunneling is partially restored. 




FIG. 3: Plot of the energies for the energy bands E-n{q) versus quasimomentum q for an optical lattice with optical depth 



II. 



OPTICAL LATTICES 



The investigations of tunneling for cold/ultracold atoms (Bose-Einstein condensates or Fermi degenerate gases) are 
based on the use of optical lattices [3 [2^. For a ID optical lattice a standing wave is created by the interference 
of two linearly polarized traveling waves counter-propagating along the a;-axis with frequency wl and wave-vector 
Al- The amplitude of the generated electric field is £{r,t) — 2fo sin(wLO When the laser detuning from 

the atomic transition is large enough to neglect the excited state spontaneous emission decay, the atom experiences 
a periodically varying conservative potential 



KiW^Fosin^ ( ^) 



(1) 
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schematically represented in Fig. [2l^a). The amplitude Vq depends on the laser detuning from the atomic transition 
and on the square of the Sq electric field amplitude [26|. The periodic potential has a spacing c^l = Al/2. This 
potential derives from the quantum mechanical interaction between atom and optical lattice photons. Therefore the 
lattice quantities are linked to the recoil momentum prcc — 27r?i/AL acquired by an atom after the absorption or the 
emission of one photon. Vq will be expressed in units of E^reo the recoil energy acquired by an atom having mass M 
following one photon exchange 

E„e = ^. (2) 

2MXl ^ ' 

Neglecting the atom-atom interactions in a Bose-Einstein condensate, our ID system is described by the following 
Hamiltonian: 

dp n f 'KX\ 

For this periodic potential the associated single-particle eigenstates in the lowest band are Bloch plane waves with 
quasimomentum q. The energies E,^(q) of the Bloch waves for the lowest bands n = 1, 2, 3 are plotted in Fig. [3]versus 
quasimomentum. Ultracold atoms are loaded into the ground state band having a minimum gap A at the edge of the 
Brillouin zone. The atomic evolution within that band or the excitation to a higher band is typically investigated. 

If a force F is applied to the atom, as schematized in Fig. ^h) , the following Hamiltonian describes the atomic 
evolution neglecting for a moment atom-atom interactions in a Bose-Einstein condensate: 



ft2 d2 



TTX 



This Hamiltonian defines the well-known Wannier-Stark problem for the electrons moving within a crystal lattice 
in the presence of an external electric field (27l - [29| . For small Stark forces F, one can picture the evolution of a 
momentum eigenstate induced by Eq. (0]) as an oscillatory motion in the ground energy band of the periodic lattice 
with Bloch period Tb [11, [H H^, where 

At stronger applied forces, a wave packet prepared in the ground band has a significant probability to tunnel at the 
band edge to the first excited band. This process of the quantum tunnel across an energy gap at an avoided crossing 
of the system's energy levels is described by the Landau-Zener tunneling [30l Isij . For a single tunneling event, the 
Landau-Zener tunneling probability is [2^ 

Plz = e I ^^.c ) , (6) 

where we introduced the Fq dimensionless force 

Fo - Fd^/Ercc (7) 

In the presence of a sequence of Landau-Zener tunneling events the Landau-Zener rate Flz to the excited band 
is obtained by multiplying Plz with the Bloch frequency — l/^B [28| . By introducing the recoil frequency 
i^icc = Ercc/h, Flz may be written 

Elz -I'rcc-Foe"^^^)' . (8) 

For the optical lattice periodic potential an alternative single-particle basis useful for describing the tunneling of 
particles among discrete lattice sites is provided by Wannier functions [iB, H^-HS HI] . The j-th Wannier function 
|j) is centered around the j lattice site, and the functions are orthonormal. In a given energy band the Hamiltonian 
for free motion on the periodic lattice is determined by hopping matrix elements, which in general connect lattice 
sites arbitrarily spaced. However, because the hopping amplitude decreases rapidly with the distance, the tunneling 
Hamiltonian may include only the J tunneling hopping between neighboring lattice sites 
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FIG. 4: Plots of the line in the space of the optical lattice depth Vo, in i?rec units, and the dimensionless force Fq dividing the 
upper (lower) region where the interwell tunneling is ten times larger (smaller) than the Landau-Zener tunneling to the upper 
band. 



where Ej defines the energy of the j-th site. For ultracold atoms in an optical lattice with depth Vq ^ -Eroci the 
nearest-neighbor tunneling energy J is given by |33i] 



In the presence of an applied force F supposing Ej = Eq = the Hamiltonian becomes 

H = FdLj2j\j >< J'l - d^' >< ^' + 11 + + 1 >< -^D ' (11) 

3 J 

However this Hamiltonian may be used to describe the atomic evolution in the ground band only when the Landau- 
Zener tunneling to the excited band can be neglected. Fig. 2] reports for a given value of the dimensionless force Fq, 
the Vo optical depth where the hopping constant J is ten times larger than Flz- 

The simulation of the temporal evolution of the Bose-Einstein condensate wavefunction is based either on the 
Gross-Pitaevskii equation based on a global mean-field description or on a many-body approach where the atomic 
number of the lattices sites is quantized [l^, [sj, |35| . Apart from the theoretical results reported in Section V, we 
will concentrate here on the mean-field approach applied to describe experimental configurations and results reviewed 
in detail in Section IV. For a realistic description of those experiments, the Gross-Pitaevskii equation was used to 
simulate the temporal evolution of the condensate wave function ip(f,t) subjected to the optical lattice and to a 
confining harmonic potential, for instance with cylindrical symmetry 

.,d , ,^ , 



2^ + Im {u^W + ^^P') + Vo sin^ (5^) +Fx + gN t)\' 



mt). (12) 



The frequencies uJx and uj^ characterize the longitudinal and transverse harmonic confinement. The atom- atom 
interactions are modeled by the nonlinear term in E g. (112I) . with the nonlinear coupling constant given hy g — 
Anfi^as/M, where Ug is the s-wave scattering length [33.l35j. Refs. [ll,[ll| introduced the g dimensionless nonlinearity 



parameter 



3 = (13) 



computed from the peak density no of the condensate initial state, to describe the nonlinear coupling relevant for 
optical lattice experiments. In the Thomas- Fermi regime of the condensate (3^. |35|. for given ojx and uJr the density 
no, and therefore g, is proportional to N^^^ where N is the number of atoms in the condensate. 
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III. RESONANT TUNNELING IN CLOSED SYSTEMS 
A. Two levels 

Quantum tunneling of a two-level system takes place in the double well potential. The quantum mechanical 
solution shows that the wave packet initially localized in one of the wells performs oscillations between the two 
classically allowed region. The period of these oscillations is related to the inverse of the energy difference between 
the symmetric and antisymmetric quantum states of the double-well system, i.e. to the energy corresponding to the 
tunneling splitting. That energy is equal to the interaction Hamiltonian between the eigenstates of the two well. 
In an asymmetric double well as that shown in Fig. [IJa) , an energy matching between a ground state on one side 
and an excited state on the other side leads to a RET between those states. The dashed lines in Fig ^a.) denote 
the eigenenergies for the symmetric and antisymmetric quantum superposition of the wavefunctions in left and right 
wells. The tunneling evolution is described by the following Hamiltonian; 

H=J2 Ei\j><j\-J{\l><2\ + \1><2\) + U (14) 

where |1 > and |2 > denote the wavefunctions of the resonant states in the left and right wells, A = Ei — E2 is the 
energy difference between the two wells, and J is the tunneling energy, U is the interatomic interaction energy and 
rij is the atom number in the left or right well. For the following analysis U represents a shift in energy of the left or 
right well. By treating at first the U = case, the atomic wavefunction may be expanded as a superposition of the 
|1, 2) states 



\^it)>= c,im, (15) 



the atomic evolution is characterized by Rabi oscillations between the two wells. For instance by supposing as initial 
condition Ci(0) — 1 and C2(0) = 0, the occupation probabilities of the left well at time t are given by 



\C,{t)f = l-\C2it)f (17) 

Therefore for the A = resonance condition of RET, a complete oscillation between the two wells at frequency 2J/h 
takes place. The atomic interaction term U shifting the -Ei=i.2 energies of the two wells may be included into the 
above equations for the occupation probabilities as a contribution to the A energy difference. Therefore the presence 
of the U interatomic energy modifies the RET condition. 

Periodic double-well structures may be created in properly chosen optical lattice or superlattice geometries. For 
cold atoms theoretical and experimental investigations were performed by a few authors [17, 36-3^. For cold atoms 
the coherence length of the atomic wavefunction is comparable to the extent of each double- well, so that the long range 
periodicity of the optical lattice plays a minor role on the tunneling properties. Therefore those investigations will be 
mentioned here. Those studies examined the new features appearing when the double-well potential depends on the 
internal atomic structure, for instance on the two electron spin states. This case was theoretically analyzed by Castin 
et al. (36l] within the context of two dimensional Sisyphus cooling. Resonant tunneling between the adjacent potential 
wells of the periodic potential for the two internal states, not present in a ID geometry, contribute with quantum 



processes to the cooling phenomena in optical lattices. Dutta et al. [37[ studied periodic well-to-well tunneling of 
^''Rb atoms on adiabatic potential surfaces of a ID optical lattice. Atoms that tunnel between neighboring wells of 
the lattice are an excellent tool for a careful study of topological potentials associated to the optical lattice. RET-like 
effects have been observed in a number of experiments to date. In Ref. fl7l|, resonant tunneling was observed for cold 
atoms trapped by an optical lattice when an applied magnetic field produced a Zeeman splitting of the energy levels. 
At certain values of the applied magnetic field, the states in the up-shifting and down-shifting energy levels were 
tuned into resonance with one another. This led to RET drastically altering the quantum dynamics of the system 
and producing a modulation of the magnetization and lifetime of the atoms trapped by the optical lattice. Hacock et 
al. [38| observed the quantum coherent dynamics of atomic spinor wave packets in the double-well potentials. With 
appropriate initial conditions the atomic system performed Rabi oscillations between the left and right localized states 
of the ground doublet, with the atomic wavepacket corresponding to a coherent superposition of these mesoscopically 
distinct quantum states. 

For ultracold atoms, Rabi oscillations in double well geometries have been investigated and measured by the 
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FIG. 5: Tunneling configuration and experimental results for the resonant tunneling of single and double atoms in a superlattice. 
The periodic double-well potential for ultracold rubidium atoms was realized by superimposing two periodic potentials with 
periodicities of Al = 765.0 nm (long lattice) and Al/2 = 382.5 nm (short lattice), and controllable intensities and relative 
phase. The depth was Vb = 12i?rcc for the short optical lattice, Vo — 9.5i5rcc for the long lattice. The upper left, lower left and 
lower right insets describe the resonant tunneling configurations for one or two atoms per well. The upper left inset describes 
the oscillating motion of the atoms between the the two wells for the conditional resonant tunneling resonance where a single 
atom can tunnel only in the presence of a second atom and the interaction energy U is matched by an applied bias. In the 
central part the amplitude of the tunneling Rabi oscillations, and the Lorentzian fit, are shown as a function of the bias energy 
A for each of the tunneling configurations represented in the insets, black dots and Lorentzian centered at A = for upper left 
one, red dots and Lorentzian centered at A = for lower left one, and blue dots and Lorentzian centered at A = 0.78-Ercc for 
lower right one. (From Foiling et al., Nature [l^. With permission by MacMillan). 

authors of [l^ [s^ . A highly parallel structure of double wells is created using optical lattice or optical superlattice 
configurations. In the superlattice configuration of fl^ the periodic potentials created by two laser standing waves at 
wavelength Al and Al/2 are applied to create a large set of individual wells. By changing the intensity of the standing 
wave lasers at the two wavelengths and their relative spatial phase, any configuration of symmetric or asymmetric 
double wells is created. In that experiment the double well investigation was performed with ultracold atoms in a 
Mott-insulator configuration having single atom occupation of the wells [l6| . The modification of the optical lattice 
potential from a periodic structure of single wells to a periodic structure of double wells, by adiabatically raising an 
energy bump within each single well, allowed to produce the asymmetric loading of each double well. 

Fig. [5] summarizes experimental results obtained in (l9| for the RET features in symmetric and asymmetric double 
wells. The tunneling of the ultracold atoms was measured as a function of the energy bias A between the wells. 
The left upper inset schematizes the case of single atom tunneling. The right lower one schematizes the tunneling of 
one atom in the presence of an energy shift produced by the atomic interaction {U term in Eq. (jl4p ). A conditional 
resonant tunneling resonance occurs, where a single atom can tunnel only in the presence of a second atom and the 
interaction energy U is matched by the bias. For these two cases the measured atomic Rabi type dynamical evolution 
between the two wells is shown the right upper inset. Because the presence of an atom in the left well shifts by U 
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FIG. 6: Unperturbed energies Ej, with j = (1.3) (in arbitrary units) of three states experiencing crossings and anticrossings 
as a function of a parameter (also in arbitrary units). Continuous lines corresponds to state |1), the dotted one to state |2) 
and the dot-dashed to state |3). The boxes marked with 1 and 3 denote regions where the tunneling is dominated by two-state 
interactions. The box marked 2 denotes a region where the three-state interaction may modify the tunneling rate between state 
|2) and |3). In region 3 without direct interaction between states |1) and |3) a locking of tunneling, corresponding to a level 
crossing with Ei = Ez, takes place. 



the level energies, a bias A = —U is applied in order to compensate the shift. Thus, a resonant tunneling condition 
is verified and the blue data denote the periodic occupation of the left well and right well, located at positions -1 and 
respectively. In the absence of an atom in the left well and without application of the bias, the tunneling is not 
resonant and the Rabi oscillations take place with a reduced amplitude and at a higher frequency, in agreement with 
the description of Eqs. and p?)) . The left lower inset schematizes the case of a correlated atomic pair tunneling, 
as produced in a second-order tunneling process. The central part of that Figure reports the amplitude of the Rabi 
oscillations versus the A bias for the different tunneling configurations, and their fits by the Lorentzian line-shapes 
predicted by Eq. (jl6p . The tunneling amplitude versus the potential bias is measured for the case of single atoms 
(black data points) and initially doubly occupied lattice sites (blue and red data points) . The blue data points and the 
Lorentzian fitted to the data point with center at A = 0.78(2)£'i.oc correspond to the conditional resonant tunneling 
resonance. The correlated pair tunneling (red circles) and the Lorentzian fit are resonant for zero bias because energies 
of both left and right wells are modified by the interaction energy U . 

While the previous description applies to single particle tunneling, quantum tunneling of macroscopic TV-body 
atomic systems introduces qualitatively new aspects to the quantum evolution of ultracold atoms, as investigated 
in [40] for Bose-Einstein condensate in a tilted multilevel double-well potential. For a double-well without tilt as 
experimentally investigated by Albiez et al. ^41.] , the so-called self trapping regimes is realized where the bosonic 
nonlinear interaction term of the Gross-Pitaevskii equation ([T^ modifies the level energies and inhibits the resonant 
tunneling between the wells. Khomeriki et al. [4^ demonstrated for a double- well structure by a pulse- wise change 
of the intermediate barrier height, it is possible to switch between the tunneling regime and the self-trapped one. 



B. Three levels 



The idea of controlling the tunneling rate between two states has led several researchers to consider the effect of 
external forces on the tunneling oscillations. Because the tunneling rate is related to the difference in the energies 
of the quantum states, a number of complicated scenarios arise when one of the states undergoes interaction with 
a third state, and that interaction may be controlled by an external parameter, for instance a magnetic or electric 
field. The tunneling wavepacket is described as a linear combination of the three initial states. Their interaction 
can drastically affect the eigenenergies of the Hamiltonian and it would be possible to explore different regimes, from 
strong suppression to enhancement of tunneling. 
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This three-level control was theoretically investigated in Refs. (43|, |44| in connection to the dynamical tunneling 
produced by time dependent potentials and for conditions as in an experiment by Raizen's group in 2001 [45[ and at 
NIST [46] . The tunneling period in the time-dependent systems is related to the differences between quasienergies of 
the Floquet states, just as the tunneling period in the time independent case has to do with the energy differences 
between the stationary states. The experimental and theoretical investigations considered the case of the tunneling 
doublet interacting with a third state associated with a chaotic region. The underl yin g classical phase space of the 
systems had a mixed regular-chaotic structure, giving the scenario of chaos-assisted f47'| or, more generally speaking, 
of dynamical tunneling [1] . 

We present here the basic of the three-level tunneling in the case of time independent potentials. Fig. |6] schematizes 
the dependence on an external parameter for the Ej energies for the \j) states, with j = 1 ... 3, in the absence 
of interactions between them. We will discuss the modifications to those energies produced by atomic interactions 
between states, supposing the presence of the interactions U12 between states |1) and |2), and U23 between states 
|2) and |3), and supposing no interaction between states |1) and |3). Notice that these interactions modify the Ej 
energies in the regions close to the energy crossings, boxes 1, 2 and 3 in the Figure, and that the tunneling frequency 
is determined by the splitting of the perturbed energies. In the box with number 1 the Ei — E2 energy separation, 
i.e. the tunneling, is dominated by the interaction between states |1) and |2). In the box denoted as 2, a three-state 
interaction takes place and the amplitude of the interaction between states |1) and |2) may be used to enhance or 
suppress the tunneling frequency between the states |2) and |3). Within the region denoted as 3, in the absence of a 
direct interaction between the |1) and |3) states a Ei — E^ crossing point exists. This crossing produces an absence 
of tunneling, this configuration being indicated as locking of the wavefunction in the initial state of preparation . 



IV. TUNNELING IN OPEN SYSTEMS 
A. Optical lattice without/with tilt 

An optical lattice is composed of an infinite number of neighboring wells uniformly distributed along one direction 
and spacing cIl = A/2 between the minima, where A is the wavelength of the standing wave laser required to create 
the periodic potential for the atoms 14] . This configuration corresponds to Fig. [2ja) . The tunneling in this system 
has strong similarities to the double- well discussed above, when the presence of physical boundaries, as in the physical 
realizations, plays no role. 

For a more general treatment we consider the case where an applied external force F produces an energy difference 
Fdi, between neighboring wells, see Fig. HJb). The atomic evolution may be studied by considering the localized 
Wannier wavefunction \ i > and the perturbations originating from the atomic occupation in neighboring sites [12]. This 
approximation is valid when the overlap of atomic wavefunctions introduces corrections to the localized atom picture, 
but they are not large enough to render the single site description irrelevant. The H Wannier-Stark Hamiltonian 
determining the atomic evolution in the absence of the interatomic interactions U is given by 

H - -^E >< + 1| + + 1 X -^D + F^lJ^jIj X j\. (18) 
i 3 

In analogy to Eq. (fTS]) the generic atomic wave function can be written as a superposition of the \ j > localized 
wavefunctions where the sum extends over all lattice sites. The temporal evolution for the Ci coefficients under the 
Hamiltonian H is given by 

ih^ = jFdLCj - J (Cj+i + , (19) 

and in the following the ground state energy Eq will be supposed to be equal to zero. The solution of these coupled 
equations with t = initial condition of atomic occupation of the i = site, i.e. Cj{t = 0) = (5j=o, leads to [i^ 



2JTr . /ttA 
sm 



(20) 



having introduced the Bessel functions Jj of j—th order. The argument of the Bessel functions in Eq. ([^0]) is an 
oscillatory function of time. Tr represents the recurrence time for the evolution of the atomic wavefunction. For 
the present case of the resonant tunneling modified by the presence of a force F, Tr = Tb whence the recurrence 
time coincides with the Bloch period Tb defined in Eq. ([5]) and is inversely proportional to the applied external force. 
The temporal recurrence of the atomic wavefunction is shown in Fig. [7] for different times expressed in units of Tr. 
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FIG. 7: Temporal recurrence of the occupation probability |Cnp versus the n position of the lattice site at different interactions 
times, . From (a) to (e) interaction times are 0, 0.2, 0.5, 0.95, 1 measured in units of Tr The occupation probabilities are 
connected by lines. Notice the reduced vertical scale at the intermediate times. The plots are obtained for the parameter 
2JTR/h = 28. 

Notice that the parameter 2JTr/^ of the Bessel function determines the range of lattice sites occupied by the periodic 
wavefunction expansion. The corresponding atomic mean-square displacement is 



TTt 



(21) 



In the limit of Fdi^ ^ J the mean-square displacement is largely decreased because of the suppression of the resonant 
tunneling, as schematized in Fig.[5][b). This suppression and the related Wannier-Stark localization of the wavefunction 
have been intensively discussed in the solid state physics theoretical literature [H, . Korsch and coworkers [13, [U 
have considered the case of an atomic distribution not initially concentrated on a single site, and instead described 
by a Gaussian distribution with root mean-square ctq- For that case the temporal evolution of the mean-square 
displacement is given by 



TTt 



cos 



2TTt 



2e 



sm 



TTt 



(22) 



In the absence of external force, taking the limit of F — >■ 0, we recover the result of a diffusion process for the 
atomic wavefunction 



2Jt 



V2Jt 



(23) 
(24) 
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B. Photon-assisted tunneling 



The above analysis can be applied also to the photon-assisted tunneling occurring when the ground states of 
adjacent potential wells tuned out of resonance by the Fdi^ static potential are coupled by photons at frequency lo as 
schematized in Fig[5fc). When the photon energy bridges the gap created by the static potential, tunneling is (partly) 
restored. The resonant tunneling is restored by a photon-assisted process when the energy provided by n photons 
matches the separation energy Fdj^ between neighboring wells. The energy resonance condition for the frequency 
is given by 

nhujn = FdL (25) 

with the integer n denoting the order of the photon-assisted resonance. This resonance may be expressed as — 
2'Ki'-Q/n in terms of the Bloch frequency. The frequency detuning from the resonance is Aw = uj — wr. 

In solid state systems, the photons are typically in the microwave frequency range and the static potential is provided 
by an electric bias field applied to the structure. Photon-assisted tunneling has been observed in superconducting 
diodes [s^] , semiconductor superlattices 0, 0| and quantum dots ^5|, |56| . 

For the photon assisted tunneling of cold and ultracold atoms, a theoretical analysis w as p erformed by Eckardt 
et al. [131 and by Kolovsky and Korsch [1^, with experiments performed by Sias et al. j59| . Ivanov et al. [60j . 
Alberti et al. (6l| and Haller et al. [63| . In these experiments a periodic time-dependent potential was applied to the 
cold atoms through a periodic spatial oscillation of the optical lattice minima/ maxima, to be referred to as shaking 
in the following. In the lattice reference frame such a backward and forward motion of the periodic potential at 
frequency uj ~ wr along one direction is equivalent to a periodic force Fi^ cos(a;i) applied to the atoms. Thus using 
the localized Wannier wavefunction introduced above for a deep lattice the atomic evolution is determined by the 
following Hamiltonian: 

i^shaking = -JY1 (I^' >< J' + 1| + \J + 1 >< j\) + [P^L + A' cos [ujt]]^ ]\] >< j\, (26) 

j j 

once again not including the U interaction term. Here K = F^di^, denoted as shaking amplitude, is the shaking energy 
difference between neighboring sites of the linear chain associated to the shaking. The theory of ref. [Irj predicts that 
when the driving takes place at the frequency 3> J/-Ercc and the resonance condition of Eq. (|25p is satisfied, the 
shaking leads to an effective tunneling rate 

Jeff(i^,C^R) = JJ„(-^). (27) 

Therefore a modification of the tunneling rate is obtained when the ratio of the rescaled shaking amplitude K = F^di^ 
and the shaking frequency times h is varied. In the experimental realization (soj the shaking frequency was fixed and 
the shaking amplitude was scanned to verify the relation of Eq. ((27|) . 

The previous analysis for the evolution of the atomic wavefunction under resonant tunneling can be applied also 
to the photon-assisted tunneling by using the approximation of a resonant dynamics introduced by Thommen et 
al. [l^ or equivalently by restricting our attention to the resonant Floquet states [G^I- In the presence of a driving 
at frequency uj and taking into account the static energy difference Fdi^ between neighboring wells, we write for the 
atomic wavefunction 

|*(t) ^.^^^^^-^{3Fd^+n^h^)t/h^■ (28) 

j,™ 

where the j index labels the well and the m index the component in the Floquet spectrum. For w close to the n-th 
order resonance condition we may restrict the terms to the resonant ones in two sums of the above expansion 

|vE'(t)>=^e-'^"^-*C'f|j>, (29) 

j 

where we have simplified the notation introducing the resonant coefficients C". 

The temporal evolution of the (7" is described by an equation similar to Eq. where Jeff determines the tunneling 
energy of the n-th order resonance. Therefore, for the photon-assisted tunneling, the occupation of the j-th lattice 
site and the mean square displacement of the atoms are the analogues to those derived previously 



m? > v^JcffTR 
dh Trh 



sm 



TTt 



2JcffTR, . / nt 



(30) 
(31) 
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with the recurrence time for this process given by 

Tr = 27r/Aa;. (32) 



This recurrence process was named as super-Bloch oscillations in refs. [58|, \6^. For the resonant case Aw = 0, the 
mean-square displacement is given by Eq. ([M)) and the occupation probabilities are given by Eq. 1^^. Notice that 
for both Wannier-Stark localization and photon-assisted tunneling, the mean-square displacement and the occupation 
probabilities have the same functional dependence if we introduce a unifying parameter for the detuning from the 
resonant tunneling. This parameter is Fdj^ for the case of an applied external force and hAuj for the case of the 
photon-assisted tunneling. Thus, the data of Fig. [7] applies also to the occupation probabilities in the photon assisted 
tunneling. 

Few experiments on optical lattices have verified or made use of the theoretical predictions of this Section. In the 
following the experiments will be characterized by the depth Vo of the optical lattice expressed in units -Ercc, and the 
photon-assisted frequency detuning Awq- 

The linear time dependence of atomic mean-square displacement predicted by Eq. (I24p in the conditions of -F = 
was applied by Lignier et at. [65l [ to measure the J tunneling energy and to verify that the experimental procedure 
reproduced the J dependence on the lattice depth Vo predicted by Eq. (fTO|) . The photon-assisted tunneling exper- 
iments [59I l60j made use of that linear dependence to measure the effective tunneling rate. In these experiments 
the linear dependence was tested for a total time larger than ten thousand tunneling times. Notice that in all these 
experimental observations the initial distribution of the atomic wavefunction was not concentrated on a single well as 
in our theoretical analysis and instead covered several wells. Nevertheless a Gaussian convolution of the initial wave- 
function spread and of the linearly expanding mean-square displacement represented a good fit of the experimental 
observations, even at earlier times where the initial width is comparable to the tunneling spread. 

The Wannier-Stark localization of the atomic cloud in the presence of an applied force F was examined by Sias 
et al. [59|] as a reduction of the mean-square displacement increasing the force amplitude at a given interrogation 
time. Fig. HJa) reports the temporal dependence of V < m? >/dh as predicted by Eq. (|2T|) . at different values of the 
parameter Fdi^/ J scanned in that experiment within the interval (0,1). In order to provide a unified description the 
time is measured in units of Tb. It appears that \J<mF> is periodic in time with period Tb while the amplitude 
of the oscillation decreases with the force until the Wannier-Stark localization regime is reached where the atomic 
motion is blocked. Fig.[SKb) shows the amplitude of the oscillation predicted by Eq. ([M)) versus the Fdi^/ J parameter. 
By comparing this dependence to the Lorentzian one occurring for a two-level system of the previous Section, it 
appears that for an infinite systems of wells the oscillation amplitude decreases more rapidly increasing Fdi^/J. For 
different values of the applied force, the maximum of the oscillation occurs at a different value of t. Therefore the 
experiment of (Hoj that measured the oscillation amplitude at a given interaction time, obtained results similar to 
those of Fig. [IJb), not precisely fitted by the inverse law as sketched in Fig. [8] 

For the photon-assisted tunneling the functional dependence on time of the wavefunction spreading on the lattice 
and the mean-square displacement was measured in l61| for a total time equivalent up to seven recurrence times in 
the case of a drive detuned by Aijj/2tt — ±5 Hz and up to one recurrence time for the Auj/2tt — ±0.260 Hz detuning. 
The measured sinusoidal evolutions are in reasonable agreement with the sinusoidal function predicted by our model 
and represented in Figs. [7] and W[a). Our model does not take into account the initial atomic distribution over several 
optical lattice sites, and in j61|, because the atomic de Broglie wavelength was shorter than the lattice period, the 
coherence degree among adjacent Wannier-Stark eigenstates was negligible. The quantum-mechanical evolution of 
the atomic wavefunction under the tunneling Hamiltonian described by our analysis is limited by the presence of de- 
coherence processes, and in (6ll | a decoherence time of 28 seconds was measured. It would be interesting to investigate 
theoretically the role of a decoherence process on the tunneling evolution. 

For the photon-assisted tunneling the mean-square amplitude dependence on the detuning Aw is given by Eq. ((2T|) 
with Tr = 27r/A(W. That functional dependence predicts that the fullwidth of the resonance line-shape Awpw, defined 
by the first zeros of the sin function, is determined by the experimental interrogation time T 

7r 

AwFw = 7^- (33) 

For interrogation times between 0.5 and 2 seconds of the experimental investigations line widths in the few Hertz range 
were measured. In the investigation of [i^l where the external force was gravity, the measurement of the resonance 
frequency for the photon-assisted tunneling with the accuracy reached by the above interrogation time allowed those 
authors to measure the gravity acceleration with ppm resolution. This shows that sensitive RET effects have a great 
potential for applications, e.g. for precision measurements. 

The recurrence process of super-Bloch oscillations was recently investigated by Haller et al. [g^l for Vb /i?roc values 
in the 3-7 range, and Au)/2tt in the 0.1=2 Hz range. The recurrence oscillations were measured up to 2.5 seconds. 
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(a) (b) 

FIG. 8; (a) mean-square displacement versus time for different values of the unified RET energy mismatch, Fdi^/J for the 
Wannier-Stark localization and h/^Lo/J for the photon assisted tunneling. Results for values 0.2, 0.4, 0.6, 0.8 and 1.0 of the 
detuning parameter, with the displacement maximum decreasing at higher values. The time dependence of \/ < m? > is periodic 
in time with period Tb. In (b) the maximum of the mean-square displacement is plotted versus Fdi^/J. The mean-square 
displacements are measured in units of the dL lattice spacing. 
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FIG. 9: Schematic of the RET process between second nearest neighboring wells, i.e. for Ai = 2. The tunneling of atoms 
is resonantly enhanced when the energy difference between lattice wells matches the separation between the energy levels in 
different potential wells. 



C. RET in optical lattices with tilt 



In spite of the fundamental RET nature and of its practical interest, for a long time the ex per imental observation 
was restricted to the motion of electrons in superlattice structures [6|. In 2007 Sias et al. [23i| observed resonant 
tunneling using Bose-Einstein condensates in accelerated optical lattice potentials. The nearly perfect control over 
the parameters of this system allowed the authors to prepare the condensates with arbitrary initial conditions and 
also to study the effects of nonlinearity and a loss of coherence. Such observation can be generalized to studying 
noise and thermal effects in resonant tunneling and underlines the usefulness of Bose-Einstein condensates in optical 
lattices as model systems for the solid state. 
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FIG. 10: Resonant tunneling in the linear regime. Shown here is the tunneling rate from the lowest energy bands of the lattice 
as a function of the normalized inverse force F^^ for Vo — 2.5 Erec lattice depth. The straight line represents the prediction 
of the Landau-Zener theory. Inset: Deviation from the Landau-Zencr prediction of Eq. Adapted from Sias et al. [23| . 

Copyright 2007 of American Physical Society. 



A schematic representation of resonantly enhanced tunneling is shown in Fig. |9l In a tilted periodic potential, 
atoms can escape by tunneling to the continuum via higher-lying levels. The tilt of the potential is proportional to the 
applied force F acting on the atoms, and the tunneling rate T^z can be calculated using the Landau-Zener formula of 
Eq. (|8]). The actual rates can dramatically deviate from Eq. (|8]) when two Wannier-Stark levels in different potentials 
wells are strongly coupled owing to the accidental degeneracy of Fig.lHlwhere the tilt-induced energy difference between 
wells i and i + Ai matches the separation between two quantized energy levels, as pointed out for cold atoms by [66| . 
Indeed, the tunneling probability can be enhanced by a large factor over the Landau-Zener prediction (see theoretical 
and experimental results of Fig. [TU)) . 

By imposing an energy resonance between the Wannier-Stark levels in different wells of an optical lattice shifted 
by the potential of the external force, one finds that the energy degeneracies occur at the values F at which FdL^i 
{Ai integer) is close to the mean band gap between two coupled bands of the F = problem 0, HI] . The actual peak 
positions are slightly shifted with respect to this simplified estimate, because the Wannier-Stark levels in the potential 
wells are only approximately defined by the averaged band gap of the F = problem, a consequence of field-induced 
level shifts 12811. 



1. Linear regime and decay rates 

Although the finite and positive scattering length of '^^Rb atoms means that the linear Hamiltonian of Eq. (|4]) is 
never exactly realized in experiments, the approximation of a non-interacting BEC is valid if the condensate density 
is maintained low. In that case, the interaction energy can be made much smaller than all the other energy scales 
of the system (recoil energy, band width, gap width) and hence it is negligible for the present analysis of RET in a 
condensate. 

Fig. [TU] shows the results of [23| for experimental investigations with low-density condensates and the nonlinearity 
parameter g less than w 1 x 10^^, defined as the limit of the linear regime. The tunneling rate Fi out of the 1-th band 
is shown as a function of Fq^. Superimposed on the overall exponential dependence of Fi/_fo on Ff^^, one clearly 
sees the resonant tunneling peaks corresponding to the various resonances Ai = 1,2,3,4. Which of the resonances 
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FIG. 11: In (a) real parts of the eigenenergies and in (b) tunneling rates for a lattice depth of V^o/Srcc = 10 and the Hamiltonian 
from Eq. Q. The eigenenergies and the tunneling rates are associated with two Wannier-Stark ladders or, equivalently, with 
two energy bands: ground state (thick black lines) and first excited state (thin red lines). The maxima of the ground-state 
tunneling rates correspond to Ai = 1,2,3, and 4. Reproduced with permission from Zenesini et al. '53]. Copyright Institute 
of Physics. 



were visible in the experiment depended on the choice of lattice parameters and the finite experimental resolution. 
The limit n = 3 for the highest band explored in [23l | was given by the maximum lattice depth achievable. 

By measuring the positions of the Ai = 1,2,3 tunneling resonances for different values of the lattice depth Vq, it 
appeared that the resonances were shifted according to the variation of the energy levels. For deep enough lattices, 
the resonance positions may be derived from a numerical simulation but can also be approximately calculated by 
making a harmonic approximation in the lattice wells, which predicts a separation of the two lowest energy levels 
{n — 1 and n = 2) of 



= 2^„c\/^. (34) 

V ^rcc 

By imposing the resonance condition AE2-1 = F'^^^d^Ai, the calculated F'*^^ resonance position results in good 
approximation with that predicted in refs (28l . \6T\ . 



2. Avoided crossings 



The accessibility of higher energy levels allowed an experimental measurement of the tunneling rates around RET 
conditions of two strongly coupled bands. The dependence of those rates on the system parameters was phrased into 
the frame of level crossing for states experiencing a loss rate. The modification of the level tunneling rate by the 
presence of a degeneracy may be described by a simple model of a two-level Hamiltonian with an energy separation 
e described by an energy crossing splitting e = and with a single level characterized by a decay rate |68l . :69] . Real 
and imaginary parts of the Hamiltonian eigenvalues are different for e 7^ 0, and two different scenarios take place 
with crossings or anticrossings of the real and imaginary part of the Hamiltonian eigenvalues. In one case, denoted as 
type-I crossing, the imaginary parts of the eigenvalues cross while the real parts anticross. In the second case, denoted 
as type-II crossing, the eigenvalues anticross while the real parts cross. The numerical simulations of ref. [24] pointed 
out that the large majority of the RET explored experimentally correspond to type-II crossings. As a consequence if a 
resonance takes place between the energy of the lower state and that of the decaying upper level, the tunneling rate of 
the lower state increases significantly. In addition the upper state experiences a resonantly stabilized tunneling (RST) 
with a decrease of its tunneling rate. Fig. [TTT a) shows theoretical predictions for type-II crossing and anticrossings 
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for the real parts of the eigenenergies associated with a RET configuration investigated experimentally as a function 
of the experimental control parameter, the Stark force determined by the Fq dimensionless parameter of Eq. ([T])- The 
associated Wannier-Stark states tunneling rates are shown in Fig. fllT b) as a function of Fq. The strong modulations 
on top of the global exponential decrease arise from RET processes originated by the energy crossings. The resonance 
eigenstates and eigenenergies for the non-interacting atoms described by Eq. Q were obtained in [24| by diagonalizing 
an open version of the Hamiltonian [H, [70l - [73j | . 

Experimental data on anticrossings in the tunneling rates are in Fig. [T^] taken from ref. [23l|. Although a direct 




FIG. 12: Anticrossing scenario of the RET rates, (a) Theoretical plot of ri,2 for Vb = 2.5i5roc with experimental points for Fi. 
(b) Theoretical plot of Fi,2 for Vb = 10i?rcc with experimental points for Adapted from Sias et al. [2^. 

observation of the discussed anticrossing scenario in two different levels for the same set of parameters was not 
possible, the experimental investigation compared the ground and excited state tunneling rates Fi and F2 with the 
theoretical predictions for two different parameter sets, as shown in Fig. 1121 This figure nicely reveals the anticrossing 
of the corresponding tunneling rates of strongly coupled levels as a function of the control parameter Fq around RET 
conditions. 



3. Nonlinearity 

This Section discusses how the experimental investigation of RET in tilted optical lattices are modified by the 
atom-atom interactions in the Bose-Einstein condensate. We focus on a parameter regime where the Stark force 
essentially dominates the dynamics of the condensate. Here the quantum tunneling between the energy bands is 
significant and most easily detected experimentally. The critical field values for which such excitations are relevant 
can be estimated by comparing, for instance, the potential energy difference between neighboring wells, Fdi^^ with the 
coupling parameters of the many-body Bose-Hubbard model, i.e. the hopping constant J and interaction constant U 

M- 

Our analysis will exclude the regime of Fq < J/ E^-cc ~ U /E^cc where a quantum chaotic system is realized [t^ - ItsI I . 
The origin of quantum chaos, i.e. of the strongly force-dependent and non-perturbative mixing of energy levels can be 
understood as a consequence of the interaction-induced lifting of the degeneracy of the multiparticle Wannier-Stark 
levels in the crossover regime from Bloch to Wannier spectra, making nearby levels strongly interact, for comparable 
magnitudes of hopping matrix elements and Stark shifts. 

For the regime of Fq ^ J jE^^c-, the effect of weak atomic interactions is just a perturbative shifting and a small 
splitting of many-body energy levels (tH . Itt} . In order to access the tunneling rates measured in the experiment of Sias 
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et al. |23j . we determine the temporal evolution of the survival probability Psur(i) for the condensate to remain in the 
energy band, in which it has been prepared initially. As proposed in ^| and applied in the experimental investigation, 
such a survival probability is best measured in momentum space, since, experimentally, the most easily measurable 
quantity is the momentum distribution of the condensate obtained from a free expansion after the evolution inside 
the lattice. Such probability decays exponentially 

PsuAt) = PsuAt = 0) exp (-R) . (35) 

In the absence of interatomic interactions in the Gross-Pitaevskii equation, e.g. for nonlinearity parameter g = 
in Eq. (1121) . the individual tunneling events occurring when the condensate crosses the band edge are independent. 
Hence Psmit) globally, i.e. fitted over many Bloch periods, has a purely exponential form, apart from the t — >■ limit 
[t^ . When the nonlinear interaction term is present, the condensate density decays with time too. As a consequence, 
the rates F are at best defined locally in time, and in the presence of RET a sharp non-exponential decay may 
occur, as discussed in [t^ HO]- Nevertheless, for short evolution times and the weak nonlinear coupling strengths g 
experimentally accessible {g defined in Eq. the global decay of the condensate is well fitted by an exponential 

law [Ullll 

PsuAt) = Psur{t = 0) exp i-Tnt) , (36) 

with rates r„ for the band n — 1 (ground band), 2 (first excited band), 3 (second excited band), in which the atoms 
are initially prepared. 

We start our study of the tunneling rate in presence of a nonlinearity by discussing the position of RET peaks. 
These peaks, whose positions for the single-particle evolution are studied in the previous part of this Section ITVCi are 
affected by the nonlinear interaction term appearing in the Gross-Pitaevskii Eq. ([T2|) for BEG. The RET resonances 
originate from an exact matching of energy levels in neighboring potential wells, and hence they are very sensitive to 
slight perturbations. A shift of the RET peaks in energy or in the position of the Stark force, predicted in [7l| for 
large value of the g parameter, is negligible for the experimental investigated nonlinearities g < 0.06, the resonance 
shift corresponding to the extremely small AFq < 5 x 10~^ value [71]. 

The g ^ 1 X 10~^ regime was entered by carrying out the acceleration experiments in radially tighter traps (radial 
frequency > 100 Hz) and hence at larger condensate densities. Fig. [T3l a) shows the Ai = 2 and Ai = 3 resonance peaks 
of the ground-state band {n — 1) for increasing values of g, starting from the linear case and going up to g « 3 x 10^^. 
As the nonlinearity increases, two effects occur. First, the overall (off-resonant) level of Fi increases linearly with g. 
This is in agreement with earlier experiments on nonlinear Landau-Zener tunneling (82l . Issj and can be modeled by a 
condensate evolution taking place within a nonlinearity-dependent effective potential Vee = Vb/(1 +4^) [111. Second, 
with increasing nonlinearity, the contrast of the RET peak is decreased and the peak eventually vanishes, as evident 
from the different on-resonance and off-resonance dependence of the tunneling rate as a function of the atom number 
N (and hence the nonlinearity), c.f. Fig. [131(b). 

The critical value of g for which the nonlinearity affects the resonance peak is estimated by comparing the width of 
the RET peaks of a band n (which essentially is determined by the tunneling width Fn^-i of the band into which the 
atoms tunnel) with the energy scale of the nonlinearity. In the experimental investigation of Sias et al. [2^ atomic 
nonlinearities corresponding to this order-of-magnitude argument were reached. For the parameters of Fig. [TU] and 
[TSf a) and the RET peak with Ai — 2, the typical width F2 of the decaying state to which the atoms tunneling energy 
is of the order of 0.2 . . .0.5 x E'rcc- Since g reflects the nonlinearity expressed in units of 8 x E'rec, this means that 
substantial deviations from the linear behavior are expected when g > 0.025 . . .0.06. The experimental observations 
confirmed that this threshold is a good estimate for the onset of the destruction of the RET peak, observed to occur 
around g = 0.02 in Fig. [T^a). 

The role of nonlinearity on the time evolution of an Wannier-Stark state localized in a single site of the optical 
lattice was also studied by Krimer et al. They predict that the nonlinearity strength leads to different regimes, 
where the nonlinearity induced shift in the energy of the lattice may enhance or inhibit RET. 

V. MANY-BODY TUNNELING 

In state-of-the-art experiments the interatomic interactions can be tuned by the transversal confinement and by 
Feshbach resonances [l^|, resulting in strong interaction- induced correlations. A good starting point for the discussion 
of true many-body effects is to use a lattice model, as introduced above for a single particle, c.f. Eq. pT|). and 
widely used in the context of strongly-correlated ultracold quantum gases [l6| . Such a lattice description has the 
great advantage that the number of degrees of freedom automatically is bounded as compared to field theoretical 
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FIG. 13: Resonant tunneling in the nonlinear regime, (a) The tunneling rates for Ai = 2 from the lowest energy band of the 
optical lattice as a function of the normalized inverse force F^^ for a lattice depth Vo/ Ercc = 3.5 and different values of the 
nonlinearity parameter, g « 0.01, 0.022, 0.033 from bottom to top. The continuous line is the theoretical prediction in the linear 
regime. The dashed lines connect the data obtained at large g values. As the nonlinearity increases, the overall tunneling rate 
increases and the resonance peak becomes less pronounced, (b) Dependence of the tunneling rate on the nonlinear parameter g 
at the position F^^ — 0.71 (solid symbols) of the RET spectrum peak and at Fq^ = 0.60 (open symbols) a the RET spectrum 
local minimum, for Vo/Eicc = 3.0. Adapted from Zenesini et al. 




FIG. 14: Sketch of most of the terms of the Hamiltonian (|38p . This model can be used to fully describe RET, since it contains 
excited levels in each potential well, in contrast to the effective model of Section FV Al 
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approaches (see, e.g. [86|,|87| and Refs. therein), and one can use it for practical numerical simulations. 

Using a single-band model, the regime of strong correlations in the Wannier-Stark system was addressed in [7J, [76 
[tsI . Issj l. revealing the sensitive dependence of the system's dynamics on the Stark force F. The single- band Bose- 
Hubbard system of [t^ [88| is defined by the following Hamiltonian with the creation aj , annihilation a/ , and number 
operators nf for the first band of a lattice with sites I = 1 . . . L: 



1 = 1 



FoEicclnf 



Ja 
2 



[al.ai + h.c.) + ^nfint - I) + taul 



(37) 



where the last term describes the on site energy. 

In order to describe interband tunneling and phenomena related to those discussed in the previous Section IV, 
such a model has to be extended to include at least the equivalent of two single-particle energy bands (as plotted in 
Fig. [3]). In the presence of strong interatomic interactions parameterized by U terms, the single-band model of Eq. 
P7p should be extended to allow for interband transitions, as e.g. realized ai F^ = Q in experiments with fermionic 
interacting atoms [21]. Doing so, the authors of [ttI. l89l] arrived at the following full model Hamiltonian for a closed 
two-band system schematically sketched in Fig. [T31 



Hit) ^ eaJ2 + £b X/ energy 
1=1 1=1 

L 

+ FgDE'rcc ^(&!a; -I- h.c.) forcc Coupling 

1=1 

- Jq ^^(c'^a|^j^a; + h.c.) + - Jf, ^^(e'^&|^j^5/ + h.c.) hopping in the bands 

/ 

1 ^ 1 
-f -J7a ^ (nf - 1) + -JJf, ^ nf (nf — 1) onsite interaction 



1=1 



1=1 

L 



1 = 1 



2Ux'^^nin\ + —Ux'^^^{b\h\aiai+\i.c.) interband interaction, 



(38) 



1=1 



where the h index and the hi, h\ creation/annihilation operators are associated to the terms of the second band. D is 
the "dipole" matrix element between the ground and excited single-particle states in a single lattice site (measured 
in 2TT/dh length units, c.f. the appendix A of [77l | for a detailed explanation of how parameters are computed from 
the physical model). 

Within this full two-band system, two dominating mechanisms promote to the second band particles starting from 
the ground band. The first one is a single-particle coupling arising from the force term 



Hi = FoDE,, 



L 
1=1 



(39) 



where the dipole matrix element D depends only on the lattice depth Vq (measured in recoil energies according to the 
definition above, c.f., Eq. ([2])). The second one is a many-body effect, describing cotunneling of two particles from 
the first band into the second band 



bjaiQi + h.c. 



1=1 



(40) 



In Eq. (|38| the tilting terms arising from the Stark force Fq have been transformed into a phase factor Q^^^^t/^B 
for the hopping terms by changing into the accelerated frame of reference f88|. This transformation nicely shows 
that the present problem is intrinsically time-dependent. Since H(t) = H(t -\- Tb) is periodic with the Bloch period 
Tb, a Floquet analysis can be used to derive the eigenbasis of the one-period evolution operator generated by H{t). 
This trick allows also the application of periodic boundary conditions, which is reasonable in order to model large 
experimental systems, typically extending over a large number of lattice sites. The Hamiltonian of Eq. (|38|) contains 
hopping terms linking nearest neighboring wells in both bands { Ja and J;,), and terms couplings different bands at 
a fixed lattice site / either by the force presence (FqD) or by interactions (Ux)- Other terms can, in principle, be 
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FIG. 15: Rate distributions for the spectrum of an open one-band Bose-Hubbard model in (a) for Fq ~ 0.47, Ja/-Ercc = 
0.22,Ua/Ercc = 0.2,!7,/£;,ec ^ 0.1 (for system size {N,L) = (7,6)) and in (b) for Fo ~ 0.17, Ja/£rec = 0.22, f/a/B.ec = 
0.2, Ux/Ercc — 0.1 ((A'^, L) — (9, 8)). In the regime where the Stark force dominates a log-normal distribution fits well the data 
(dotted in (a)), whilst a power-law P(r) oc F"^ distribution is found with 2: ~ 2 in the strongly coupled case (dashed line in 
(b)). 



included, yet they turn out to be exponentially suppressed for sufficiently deep lattices which are well described by 
Bose-Hublaard like models p^ . 

Because of its complex form and the large number of participating many-particle states, the above Hamiltonian 
is hard to interprete and to treat even numerically, for reasonable numbers of atoms TV and lattice sites L. Two 
approximate treatments will be presented in the following. Section IV Al uses an effective one band model which 
nevertheless takes the coupling terms between the bands of (pS)) into account. While this model is valid for small 
interband couplings. Section IV Bl presents analytical and numerical results for the full model ((38)) . which on the other 
hand is valid for arbitrary interband couplings but is perturbative in the atom-atom interaction terms Ua,Ui,, and U^- 



A. Open one-band model 



Instead of using a numerically hardly tractable complete many-band model, we introduce here a perturbative decay 
of the many-particle modes in the ground band to a second energy band. This novel app roach when applied to the 
Landau-Zener-like tunneling between the first and the second band [H, [2^, [6^ [63, lU, |8^] predicts the expected 
tunneling rates and their statistical distributions. 

To justify this perturbative approach, it is crucial to realize that the terms of Eqs. ([39|) and (|40| must be small 
compared with the band gap AE = et — ea and indeed FqD, Ux ^ AS for the parameters of Fig. [151 As exercised 
in detail by Tomadin et al. [zl, [77[, from these two coupling terms by using Fermi's golden rule one can compute 
analytically the corresponding tunneling rates Ti{s) and T2{s) for each basis state labeled by s. Those rates allow the 
computation of the total width F(s) = Ti{s) + T2{s) defined by the two analyzed coupling processes for each basis 
state \s) of the single-band problem given in Eq. p7p . The T{s) are inserted as complex potentials in the diagonal 
of the single-band Hamiltonian matrix. Alon g w ith the statistics of the level spacings defined by the real parts of its 
eigenspectrum Re{Ej} studied in (t^. [76l - l78i l88l| . the statistical distributions of the tunneling rates Tj = —21m Ifil 
may be analyzed, as done in Fig. [T5l For the regime where the motion of the atoms is localized along the lattice [2^ 
that distribution is in good agreement with the expected log-normal distribution of tunneling rates (or of the similarly 
behaving conductance) 19011. I n that regime the Stark force dominates and the system shows nearly perfect single- 
particle Bloch oscillations , the distributions agreing with those predicted from the localization theory [q^, [sH . 
On the other hand, when the Stark force is comparable with Ja and Ua and all modes of our Bose-Hubbard model are 
strongly coupled, the rate distribution of Fig. llSf b) follows the expected power-law for open quan tum chaotic systems 
in the diffusive regime [9l| . This regime shows strong signatures of quantum chaos [TJ, ItgI - ITSI [ssj , which manifest 
also in the rate distributions [zl, [T^I • 
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FIG. 16: Population in the upper band as a function of time for the rescahng parameter a = (black dotted line), 0.2 
(faint red line) and 0.5 (green thick line) in a closed two-band model. Clearly visible are the interaction induced coUapes 
and revivals of the RET oscillations between the bands. Other parameters are _Fo = 1.87 (dominating energy scale!) and 
Ja = 0.1, Jb = 0.77, Ua = 0.023, Ub = 0.014, = 0.01, tb - ta = 3.38 (aU in recoil energy units) and D = -0.16 in length units, 
and {N,L) = (5,6). 



B. Closed two- band model 



Since the model introduced in the previous Section IV Al cannot account for resonant tunneling between a ground 
level of one well and an excited level of another well, a different model which applies also for strong transitions between 
the bands was investigated by Plotz et al. [s^. This model is based on the full Hamiltonian of Eq. (155)) sketched 
schematically in Fig. [Ml 

When the Stark force is tuned to the value where RET occurs for the single particle problem (c.f. Section flVCp . 
the strong coupling of the atoms prepared in the ground band into the excited band plays an important role. Since 
the model is closed, i.e. higher bands are neglected, there is no asymptotic tunneling as in the experimental situation 
described in Section [IV CI As a consequence, we observe an oscillation of the probability of occupying the lower and 
upper band, respectively, which is particularly pronounced at RET conditions. For a single particle in our lattice 
model, such RET oscillations can be understood easily, since in Floquet space (remembering that our Hamiltonian 
of Eq. ([551) is periodically time-dependent) the problem reduces to an effective two state model of resonantly coupled 
states [8^, [H^I • In this effective description, the evolution corresponds to the two level Rabi problem of quantum optics 
[9^ . For non- vanishing atom-atom interaction, the situation complicates, of course, and we expect a degradation of 
those single-particle Rabi oscillations. This is illustrated in Fig. [121 The period of the single-particle interband 
oscillation is given by the following formula derived in [sl] : 

iosc 1 ^^^^ 



ch 



2DJA^ 



Jb- 



where Ai is the resonance order introduced in Section [IV CI and Jai the Bessel function of the same order. 
For a Stark force Fq not satisfying the RET conditions, the coupling to the upper band is strongly suppressed, and 
almost negligible at least for small particle-particle interband interactions Ux- On the other hand if Ux dominates, 
strong interband coupling is possible even for small forces Fq. The latter strongly correlated regime of two energy 
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FIG. 17: CoUape and revival times extracted from data (symbols for two different system parameter sets) as shown in Fig. [16] 
versus the inverse of the atom-atom interaction rescaling factor a. As expected for a two-state Rabi problem perturbed by a 
coupling to additional states, both times scale inverse proportionally to a. The dashed lines should guide the eye. 



bands is extremely hard to deal with, especially if one is searching for analytical predictions for the interband dynamics. 
The results shown in Fig. [12] are just a small step in this direction. In the limit of small atom-atom interactions, 
the observed collapse and revival times can be determined analytically in good approximation. We quantify small 
interactions by artificially rescaling the parameters Ua,Ub,Ux, which would be obtained by a given scattering length 
and a given depth of the optical lattice potential by a constant factor < a < 1. For the results of Fig. [TG] a 
was chosen to be zero (black dotted line), 0.2 (faint red line) and 0.5 (green thick line). The analogy with the Rabi 
oscillation problem even carries over to those values of interaction strength, since we observe a collapse and later on 
a revival of the periodic oscillation of the population. Collapse and revival time scale inversely proportional with the 
strength factor a, as shown in Fig. 117] where the revival time is well approximated by the formula derived in [soj 

^revival 2_fo fAO\ 

^ " aU..P, (A) (A) ' 

with the zcroth order Bessel function Jq. This formula arises from a perturbative calculation of the effect of atom- 
atom interactions for small aUa.b,x ^ Fq starting from the single-particle solution, which itself is known within the 
effective two-state model, and assuming a delocalized initial state along the lattice. From Eq. (|42]) the collapse time 
was estimated in j9^] as icoiiapsc ~ trcvivai/ i'^o's), with the effective number CTs of additionally coupled many-particle 
states as compared to the single-particle two-state model. This collapse is analogous to that of the Rabi oscillations 
in the presence of atomic interactions, or to the collapse arising whenever the phase evolution of each s basis state is 
nonlinear in the particle number. Notice that the collapse and revival phenomena of Fig. [16] stem from a degradation 
(arising from interactions) of single particle interband oscillations (with original period given by Eq. (|4ip which jus t 
depends on the force Fq). So, even if there are analogies to the collapses and revivals observed in BEC [isl. l93 - l96l | . 
their origins are different. In the BEC investigations the collapse-revival oscillations were produced by the interaction 
within a single-band (in [96| by atomic interactions depending on higher power of the well occupation number). 
Therefore those oscillations would not at all occur when the lower band nonlinear interaction(s) is (are) suppressed, 
equivalent to Ua = in the model here discussed. 

The above steps may be expanded in different directions within the realm of true many-body dynamics and 
tunneling, with great perspectives for many-body induced RET effects. Remaining questions are, for instance, the 
study of the strongly correlated regime of strong particle and strong interband interactions simultaneously, and the 
enlargement of our closed two band model in order to allow for a realistic description of experiments similar to the 
ones reported in Section [IV CI now carried over into the realm of strong many-body interactions. 
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VI. CONCLUSIONS AND PERSPECTIVES ON RET 

This chapter has presented and discussed the RET investigations performed with cold and uhracold atoms. Owing 
to the reached high level of control on the atom initial preparation and on the realization of potentials with arbitrary 
shapes, the atomic physics community has reproduced and analyzed basic quantum mechanics phenomena well estab- 
lished, and with important applications, within the solid state physics community. An important feature associated to 
the investigations on the atoms, compared to those on electrons in a solid, is the absence of decoherence phenomena. 
Therefore quantum interference phenomena may play an enormous role on the tunneling temporal evolution of the 
cold atoms. For the ultracold atoms an additional characteristic is the presence of interatomic interactions, that 
modify the position of the energy levels and therefore greatly influence the RET. In more complex configurations 
the atomic interactions lead to a very complex Hamiltonian whose action on the atoms requires large computational 
efforts or analyses based on perturbation approaches. 

Our analysis was restricted to potentials which are either not explicitly time-dependent or lead to a temporal evo- 
lution of the atomic wavefunction corresponding to an adiabatic evolution of the atomic system. Tunneling processes 
produced by a non-adiabatic atomic evolution are described in other chapters of this volume. 

Macroscopic quantum tunneling is an important direction of research well investigated by the solid state physics 
community. Up to now no clear evidence of that tunneling was reported by the BEG community even if configurations 
for the occurrence of macrosco pic quantum tunneling in Bose-Einstein condensates have been proposed by different 
authors. Ueda and Leggett (97l.l98|| examined the instability of a collective mode in a BEC with attractive interaction 
induced by macroscopic tunneling. Thus, a collective variable the spatial width of BEC is analyzed a as a tunneling 
variable. Carr et al. [80| studied BEC in a potential of finite depth, harmonic for small radii and decaying as a 
Gaussian for large radii, which supports both bound and quasi-bound states. The atomic nonlinearity transforming 
bound states into quasi-bound ones, leads to macroscopic quantum tunneling. The experimental observation of such 
macroscopic tunneling would enlarge the quantum simulation configurations explored with ultracold atoms. 
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